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Abstract 

This document reviews the general approach to correcting the process e + e~ — > X° — > // for 
radiative effects, where X° represents an exchanged boson arising from some new physics. The 
validity of current methods is discussed in the context of the differential cross section. To this 
end the universality of the dominant QED radiative corrections to such processes is discussed 
and an attempt is made to quantify it. The paper aims to justify, as much as possible, the 
general approach taken by e + e~ collider experiments to the issue of how to treat the dominant 
radiative corrections in fitting models of new physics using inclusive and exclusive cross section 
measurements. We conclude that in all but the most pathological of new physics models the 
dominant radiative corrections (QED) to the tree level processes of the standard model can 
be expected to hold well for new physics. This argument follows from the fact that the phase 
space of indirect new physics searches is generally restrictive (high events) in such a way 
that the factorization of radiative corrections is expected to hold well and generally universal 
infrared corrections should be prevalent. 
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1 Introduction. 



In this document we discuss how radiative corrections which significantly affect the standard 
model difermion process maybe exploited such that they can be used with analogous difermion 
processes containing new physics components, i.e. we discuss the universality of radiative 
corrections to the difermion process. To this end we must derive again the significant QED 
corrections to the difermion process. The results of our derivations agree with those in the 
literature. It has also long been known that QED radiative corrections can exhibit soft behavior 
which is independent of the underlying process, the so called hard process (see e.g. [HE]), it 
is this behaviour that we wish to find and study. The known radiative corrections generally 
manifest themselves as long and complicated analytic expressions and one generally has no clue 
as to whether they are significant or universal or, ideally, both. In reviewing and "dismantling" 
known results we are able to thoroughly check the dependence of radiative corrections on the 
hard scattering process thus highlighting, as much as possible, universal parts of radiative 
corrections which are relevant to more general difermion processes. 

2 Fitting Beyond the Standard Model Physics. 

Currently there are a large number of models involving physics beyond the standard model 
which may be constrained from measurements of the difermion process e + e~ — > ff. New physics 
models appear frequently in the literature which propose modifications to these observables. 
It is common for these models to involve the exchange of a new gauge boson mimicking the 
exchange of a Z° or photon in difermion processes. Let us generically denote these particles 
X°. The result of such exchanges is to increase the measured cross sections above that which 
one would expect for difermion processes as defined above, the differential cross section will 
also be modified. 

The new physics contribution to the measured difermion cross section is always given as 
a function of some parameters of the theory, typically the energy scale of the new physics. 
The total new physics contribution due to exchange of X° bosons comprises of a term due to 
interference of X° exchange amplitudes with each of the Z° and 7 exchange amplitudes, as well 
as a pure X° exchange term. Given that new physics effects are weak the largest contribution 
from them comes from the interference term. By adding this new physics prediction to the 
prediction for genuine difermion events and fitting the result to the data for various values of 
the new physics parameters one can obtain x 2 as a function of the parameters. From this x 2 
function one can obtain confidence limits on the parameters. 

Due to the fact that the standard model generally agrees well with experimental measure- 
ments any new physics effect which may be occurring in current e + e~ collider experiments must 
be small. Generally one is trying to set the most stringent constraints on the new physics model 
one can. It is for these reasons that when selecting events with which to set limits one should 
not include events for which <C \/s^ as previous runs of the experiment will have directly 
probed these energies (and found nothing) so including such effects will essentially dilute the 
effects of a new physics contribution resulting in a more relaxed limit. Obviously there is a 
trade off between the number of events in the sample and the amount by which new physics 
would be diluted. At LEP2 limits are set on the various models of new physics using only 
samples of events for which ^fs' > 0.85a/s. 

Given that radiative corrections have striking consequences for difermion events it would 
be naive to think that they are inconsequential for new physics processes. New physics cross 

Wyfp is the so-called reduced invariant mass. Vs^ < s/s due to bremsstrahlung emitted from the initial state 
particles. Generally ^/s' is the invariant mass of the exchanged boson, it is defined in appendix|Hl 
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sections and differential cross sections that appear in the literature are only ever given to tree 
level, to obtain realistic limits for the parameters of the model one should try to correct these 
predictions for radiative effects. Effective theories are generally non-renormalizable, so the 
idea of computing radiative corrections to processes derived from them may not seem sensible. 
Nevertheless the idea of fitting the born level predictions is certainly not correct. Finally, in 
the event that explicit calculations of radiative corrections were possible, to do this would be 
unfeasible given the current rate at which new physics scenarios are conceived. 

The conventional method employed by experiments to improve the tree level new physics 
predictions of an observable O^p, to include radiative effects 0'\ NP , is to multiply the tree level 
prediction by a simple factor of the standard model prediction including radiative corrections 
0'\ SM divided by the standard model prediction at tree level 0\ SM , 

0'\ NP =0\ NP x^. (2.1) 

This correction factor embodies the aggregated effects of all of the standard model radiative 
corrections. This approach amounts to claiming that all standard model radiative corrections 
factorize trivially from the born level quantity, that radiative corrections are independent of 
the born level process. In the rest of the paper we investigate whether the approach of equation 
12.11 is at all valid. 



3 Radiator Functions: A "Black Box" Empirical Study. 

The effect of bremsstrahlung on e + e" collisions is highly significant, the data show that the 
cross-sections for events with yfs! > 75 GeV is larger than those events with > 0.85 by a 
factor of two or more. If one imagines a single bremsstrahlung correction to a difermion process, 
in which one of the colliding particles emits a photon, one effect will be to lower the invariant 
mass of the basic difermion reaction we are trying to measure from y/s to ^fs' (or \ft to \ft' 
for t-channel processes). Assuming that the processes of emitting a photon and undergoing a 
difermion reaction are otherwise independent, that is to say they factorize, one can think of 
the cross section a (s)measured at some value of s/s, in some range s m j„ < s' < s as being 
comprised of a weighted sum over all of the tree level cross sections axree (s') in that y/s' range, 

a(s)= I ds' R(s,s')a T ree(s') . (3.1) 

" $min 

R (s, s') is the weighting function, it is known in the literature as the radiator function. Given 
that the cross section is proportional to the number of events in the sample we could write 

N (s) = f ds' R (s, s') N Tree (s') , (3.2) 

$min 

which tells us that ds' R (s, s') Nrree (s') is the number of tree level difermion events in the 
range s' — > s' + ds' entering the measured sample of N (s). Therefore ds' R(s,s') is the 
fraction of the tree level cross section CTree (s') that goes into making er(s). What is being 
described here is essentially the structure function approach, this is the idea that the electron 
and positron are not to be thought of as fundamental but rather they should be thought of as 
objects containing truly fundamental electrons, positrons and photons. With this in mind one 
can define a structure function for the electron D e (x) . D e (x) is a probability density function 
of x which is the fraction of the longitudinal component of the electrons momentum that it has 
when it undergoes the tree level hard scattering process. D e (x) dx is therefore the probability 
that the electron collides with the positron with a fraction x of its original longitudinal four 
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momentum. The positron structure function is the same as the electron structure function. 
Neglecting transverse momenta for the time being we have 



s' = X1X2S. (3.3) 

With these definitions the cross section for an electron and positron to undergo a difermion 
interaction such that s' > s' „•„ is 



f dxx [ dx 2 D e (x l )D e (x 2 )a Tree (x l x 2 s) . (3.4) 



Defining 1 — x — x\x 2 and demanding it be greater than imML this can be written 



a (s) = Jo s dx [Jl_ x dx 2 jr 2 D e {±=2) D e (x 2 )\ a Tree ((1 - x) s) 

s min 

= Jo s dx H (x,s)a Tree ((l - x) s) 



(3.5) 



allowing an identification between the radiator function and the structure functions. Clearly 
H (x, s) and R (s, s') are related by a factor. 

The most general new physics distribution to fit is the differential cross section. In this case 
the quantity 0\ NP to be corrected via equation 12. II is the value of the differential cross section 
in a given angular bin. Conventionally one would then use a number of correction factors 
determined using the aforementioned standard model cross sections in that bin to correct the 
bins of the new physics tree level angular distribution. It is however hard to attach a physical 
meaning to these correction factors, they do not lend themselves to a description in terms of a 
radiator function or structure functions. Ideally one would like to employ a radiator function 
approach to correcting new physics observables and so correct them by folding the given tree 
level quantities with the radiator function. This raises two questions, does a radiator function 
approach exist for the differential cross section and if it does exist how should it be implemented? 

Radiator functions for the differential cross section do exist and are present inside e + e~ — > 
f f simulation programs. There are a variety of such programs which are freely available. There 
are two kinds of e + e~ — > f f simulation packages. Semi-analytic simulations such as TOPAZO [TH] 
and ZFITTER [H] produce predictions of observables like cross-sections and differential cross sec- 
tions. They contain analytic expressions for these quantities corrected for a myriad of radiative 
corrections all in the context of the standard model. Some outputs of these programs involve 
numerical integrations hence "Semi-analytical". In addition to the semi-analytic packages there 
are difermion Monte Carlo programs such as KKfli [TJj which simulate actual difermion events 
i.e. there output includes four vectors for the various final state particles. At the heart of 
the Monte Carlo is a random number generator, as the number of generated events tends to 
infinity the cross sections etc that are calculable with the generated events tend to those of the 
standard model i.e. random fluctuations die away. Naturally the semi-analytic packages have 
the advantage of requiring much less computational time than the Monte Carlo programs and 
the ZFITTER package appears to be the most popular of these in the experimental community. 

In [HI (and the ZFITTER manual [12]) the differential cross section for difermion production 
in the presence of ISR is given in the form^ 3 ) 

da (s, cos 6) f s , , „ , , ... . s 

, V Q = ds' R(s,s',cos9)a Tree (s'). 3.6 
dcost> Js' . 



( 3 'The form of emmtion 13.61 in [Tfi| differs slightly from what we quote which is a simplified version. The 
only difference of note is that split the born level cross section into a part which is symmetric in cos 9 and 
an asymmetric part, each of which is convoluted with a slightly different radiator function. Equation 13.61 is 
appropriate for the time being. 
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Note crucially that all angular dependence is contained within the radiator function. Let us 
denote the cross section for a new physics process incorporating radiative corrections by (Jnp 
and the corresponding tree level cross section cr TreeiNP . Assuming that it is possible to extract 
R (s, s', cos 9) from the ZFITTER package 

da NP (s, cos 9) f s , , , 

= / dS R{S,S , COS 9) (J T ree,NP [S ) , (3.7) 

dcos6> '-' 



min 



will not be the correct differential cross section for a new physics process, obviously there is 
absolutely no reference to the angular distribution of new physics in the above formula, only to 
that of the standard model. Studying references |3| and [HI] in a more detail one finds that in 
fact the radiator function which exactly factorizes at the level of the integrated cross section (by 
integrated we mean integrated over the full solid angle), as shown in !3.6[ also largely factorizes 
from the differential cross section i. e. 13.61 is of the form 

= IL, M TIT " + J,l„ ds' fl 2 (.y , cos 6) a Tree (,') 

= ff ds' fli (s, s') ^T„Avo,a> + Non . factorizab i e . v ' ' 

min ^ COS o 

In 13.81 R (s, s', cos 9) has been split into a factorizable part^ 4 ^ R\ (s,s') ^ms"^ and a non- 
factorizable part R 2 (s, s', cos9). If the parts of the above which do not factorize into the 
differential cross section are negligible then it appears that the radiator function does not 
depend on the tree level differential cross section, in fact it does not depend on angles at all. If 
this is true then Ri (s, s') would appear to be independent of the tree level physics analogous 
to a structure function (the electron structure function is discussed in the next section), then 
one would appear to have a recipe for making new physics angular distribution corrected for 
radiative effects 

da NP (s,cos9) f s ,,5, da T ,Np(s,cos9) 

ds Ri{s,s) . (3.9) 



dcos# Js' - dcos# 

mm 

Assuming it is true that the non-factorizable parts of R (s, s', cos 9) are negligible we then wish 
to obtain R\ (s, s'). 

The ZFITTER package contains just this information though as much as the function is 
convoluted with the standard model in the above equations it is convoluted to an greater 
extent with the standard model in terms of the actual coding of the package. A different, more 
practical means of deconvolving the radiator function is required. Denoting a (s, q,) as the cross 
section for the difermion process occurring such that the final state fermion goes into the cos 9 
bin in the range cos 9b < cos 9 < cos 9^+1 denoted by q, we have 

a (s, c b ) = f s s , . ds' {R, (s, s') a T (s, c b ) + J^ +1 d cos 9 J s s , ds' R 2 (a, s', cos 9) a T (s') } 

min \ & mzn J 

= Jg/ ds' \ Ri (s, s') ax (s, Cb) + J c c ° s s e 6+1 d cos 9 Non-factorizable [ 

where here "Non-factorizable" is the same as in 13.81 The method by which we extract the 
radiator function assumes that the non-factorizable part is zero. If this assumption is true one 
should get the same function R (s, s') regardless of what binning one uses in the equation above. 
The fact that the non-factorizable part has a complex angular form (this can be checked by 
reference to |16| ) will mean that the R\ (s, s') which one extracts from the program could look 
very different depending on which angular bin is used and, importantly, the size of the non- 
factorizable contributions. If the non-factorizable component really is small then the Ri (s, s') 
which is extracted from the program should not be a function of the bin used to extract it, 

( 4 ) Henceforth we use the word factorizable to mean that the radiative corrections can be represented by a 
radiator function convoluted with the differential cross section da. 
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Radiator Function 



Figure 1: The radiator function Ri(s,s') with s = (200 GeV) 2 , obtained from the ZFITTER 
package by extracting it separately for each angular bin for the final state. 

it should be flat when plotted against the bin used to extract it. The author has devised a 
means of extracting a numerical version of R\ (s, s') from the ZFITTER package by making finite 
approximations to the integral over s . The resulting numerical radiator function is shown in 
figured it is, as hoped, flat across the angular region shown. 

4 Structure Functions, Leading Log Approximation. 

The purpose of this section is to add weight to the empirical observations of the last section 
by deriving the electron structure function. The idea for a structure function for the electron 
originates from the paper of Gribov and Lipatov [TT]. In the papers of Kuraev, Fadin, Jadach 
and Skrzypek [H E] the authors consider the process of e + e~ annihilation as being a Drell-Yan 
process and exploit the perturbative nature of QED to calculate the electron structure func- 
tion analytically using the Gribov-Lipatov Altarelli-Parisi equations. This structure function 
approach represents a leading logarithmic approximation (LLA), this corresponds to assuming 
near collinear emission of photons from the incoming e + e~. The radiator function in the litera- 
ture is by no means unique. The structure function approach / the leading log approximation, 
is one way to take into account the effects of ISR. The structure function approach is intuitive, 
simple and it produces a cross section to better than 0.015% this alone is a good reason 

to discuss it. For clarity and completeness we review the results of Kuraev and Fadin (KF) as 
well as Jadach and Skrzypek (JS) with the help of a combination of papers El GDI H2]- 
In doing this we explicitly confirm (as noted in the introduction of jjj) that no element of the 
underlying hard scattering process enters this description of ISR/ 6 ) and check whether or not it 
may be applied to the differential cross section as well as the integrated cross section. 

The non-singlet structure function D NS (x, s) corresponds to the probability of finding an 

( 5 )This number is quoted in relation to the accuracy of the LLA around the Z° peak at LEP1 however in LEP2 
studies the LEP1 and higher order LEP2 structure functions have been shown to produce cross sections 
agreeing to better than 0.1% for cuts like those used in obtaining our non-radiative class of events. 

( 6 'lt had been a cause of concern that, contrary to QCD, the perturbative, calculable nature of this structure 
function may somehow allow for some specifically process dependent corrections to enter it, as is suggested by 
the explicit angular dependence in the radiator functions of [Tfi| . 
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electron with longitudinal momentum fraction x 'inside' the original on-shell electron at energy 
scale \/s. The generic starting point for determining the structure function is Lipatov's equation 
(see e.g. 0) 

c» M^(i- I )+r^t i /'- f> (*o D " s (-• s • c- 1 ) 

i m 2 Si 27T Jx X\ \X\ J 



where P (z) is the regularized splitting function 

i 



P (xO = -6(1- Xl ) dz ±±^-. (4.2) 

1 — Xi io 1 — 2; 



For our purposes a more useful form of the regularized splitting function involves introducing 
a cut off at z = 1 — e where e <§C 1 is very small 

1 , 1 +z 2 „ /3 . , \ „ . 1 2 



dz = lim2 - + loge -2e + -e J . (4.3) 

o 1 - z e-o V4 & / 2 v 7 

In this way we can define a different regularized splitting function 

P(xi)= 2<y(l-xi)(| + loge)+0(l-e-xi)i^. (4.4) 

which is finite when integrated over in the limit e — ► 0. Clearly this form is the natural extension 
of !4.2[ in the limit e — > 14.21 and 14.31 are identical. Henceforth we omit writing lim e _».o an d drop 
all terms from our calculations which are vanishing in the limit e -— 0. 

Lipatov's equation 14. II has an obvious iterative solution, to see this it is convenient to first 
rewrite the second integral using 

f 1 dx,dx 2 5 (x - Xl x 2 ) P ( Xl ) D NS (x 2 , 8l ) =f^P (x0 D NS (-, 8l ) . (4.5) 

J0 Jx Xi \Xi J 

In addition we approximate the running coupling constant to a constant a (s) = a and define 
r) = 7] (s) = ^log ^2 , r)i = rj (s^ therefore 

*» = (4.6) 

Si IT 

With these substitutions Lipatov's equation 14.11 becomes. 

D NS (x, ri) = 6(l-x) + \ [ V dm f 1 dx 4 dx 2 5 (x - x x x 2 ) P (n) D NS (x 2 , m) (4.7) 

4 Jo Jo 

Substituting relation l4~6l into Lipatov's equation l4.1l and iterating it i.e. substituting D (x iy r]i)&s 
defined hv 14.11 back into 14.11 we have an expansion in 77; 



D NS (x, 77) = 6{l-x) 

+ 1 Jo d? 7i Jo 1 dx^x 2 S(x- xix 2 ) P (xx) 6(1- x 2 ) 

+ \ Jo dr)i Jo 1 dxtdx 2 S (x - x x x 2 ) P (xi) \ j^ 1 dr] 2 Jq 1 dx 3 dx 4 5 (x 2 - x 3 x 4 ) P (x 3 ) ... 
+ Oirf) 

(4.8) 

The fact that the integrals Jq i drji are nested, that is to say the upper limit of one integrand 
is the dummy variable of the one preceding it, results in a sequence of exponential factors ^ 
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multiplying the terms above. 14.71 is condensed by simple delta function manipulations 

D NS (x,rj) 
= S(l-x) 

+ lx(2)V(x) 

+ l x (l) 2 Iodx 1 dx 2 5(x-x 1 x 2 )P(x 1 )P(x 2 ) (4 ' 9) 

+ g x (4) So dxxdx 2 5(x- xix 2 ) P (xi) /J dx 3 dx 4 5 (x 2 - x 3 x 4 ) P (x 3 ) P (x 4 ) 
+ 0(r/ 4 ) 

Denoting the convolution integrals above by 

•1 

dxidxj 5(x- XiXj) f (xi) g (xj) = f (.) <g> g (.) (x) , (4.10) 



00 -r 

D NS (x,s) = 5(l-x)+J2 1 ( l i) (4-H) 



equation 14.81 becomes 

o 

n NS (t, = ^1-^ + ^ 

where 

P® 3 (x) = / dxidx 2 5 (x — x\x 2 ) P (xi) / dx 3 dx4 5 (x 2 — £3X4) P (23) P (2:4) (4-12) 

Jo 

e£c. 

The 0(?7) term is simply the regularized splitting function, ^rjP(x). In calculating the 
O (jf) term first we calculate the theta term, we proceed essentially in an identical way to that 
in which we calculated the regularized splitting function. We calculate a regularized version of 
this term by taking x to be less than a cut off x < 1 — e. 

9 (1 — x — e) Jq 1 dxidx 2 5 (x — x\x 2 ) P (x\) P (x 2 ) 
= 6 (1 — x — e) 4 (| + log e) Jq 1 dxidx 2 5 (x — Xix 2 ) 5 (1 — x x ) 5(1 — x 2 ) 
+ (1 _ x _ e ) 2 (I + log e) J*" 6 dx! Jo 1 dx 2 5 (x - xxx 2 ) 5 (1 - x 2 ) g| (4.13) 
(1 - x - e) 2 (I + log e) J 1 dx x Jq 1 " 6 dx 2 5 (x - x x x 2 ) 5(1-^)^ 



-XI 

_2 1 1 ™2 



+ 6> (1 - x - e) Jo dxi Jo dx 2 5 (x - Xix 2 ) (1 - e - a^) (1 - e - x 2 ) r^lS^ 

The first term is zero due to the cut off we impose on x, the other terms are also straightforward 
to evaluate. As stated earlier small terms O (e) and above, vanishing in the limit e — > 0, are 
dropped. The total theta term obtained with the cut off x < 1 — e is 

9(1- e-x) Jq dxidx 2 5 (x - x x x 2 ) P (xi) P (x 2 ) 
= 0(1 -x-e)2(^f (21og(l -x) -logx + §) + ± (1 + x) logx - 1 + x) ^ ' 

The delta term, analogous to the delta function term in the regularized splitting function I4.4| 
"soaks up" the divergent part of the term above arising from the limit x — > 1 (e — ► 0) 



- lim^o Jo e dx 2 (2 log (1 - x) - logx + §) + \ (1 + x) logx - 1 + x 



+ loge -C(2) 



(4.15) 



Finally the full regularized (9 (?y 2 ) term is 

I (4) So dxidx 2 5 (x - xix 2 ) P (xi) P (x 2 ) 
= (2) 2 5 (i_ x )('2(| + log e ) 2 -2C(2)) . (4.16) 



+ (2) 0(l-x- e )(i±J (2 log (1 -x) -logx + I) + ±(l + x)logx-l+x) 



The third order term, O (rj 3 ), is simply 

£ (z) Jo dxidx2 ^ ( x - XlX2 ) _^ dx 3 dx 4 5 (x 2 - X3X4) P (x 3 ) P (x 4 ) . (4.17) 

Clearly the Jq dx 3 dx 4 part of this term is given by simply substituting in the O (r/ 2 ) term with 
the replacement, x — > x 2 . Noting this simplification the resulting calculation is significantly 
simplified and a long but straightforward sequence of integrations gives the O (rj 3 ) theta and 
delta tenW 7 ) of 

The Lipatov equation can also be solved analytically in the soft limit x — > by means of a 
Mellin transformation. The Mellin transform is defined as, 

jjNS(z) / ) = f l dx X z-1 D NS , n / 41g x 
•/ 



for which the inverse is, 



I rc+ioo _ 

P^ 5 (2, 77) = — / d^ x- 2 P (2) (77) . (4.19) 

Jc— zoo 

is bounded for some z > in which case the 



The transform exists if Jq dx x z ~ l D NS (x, rf) 
inverse exists for c > z. 

Taking the Mellin transform of Lipatov's equation (14.71) gives 



4 jo 

Differentiating the above with respect to 77 gives, 

Ad m « (77) = IpWD^W (t/) . (4.21) 

Integrating this equation one finds a simple expression for the Mellin moments of the structure 
function in terms of those of the splitting function 

f)NS(z) ^ = exp Qr?P (z) ) (4.22) 

where we have used as an initial condition D NS ^ {m 2 e ) = 1. The Mellin moments for the 
splitting function are 

pW = 2( i + log £ )-^ dl __ + ____ ,4.23) 

Using the Taylor expansion = Y^= z xU the integral in !4.23l is trivial, as before we drop terms 
vanishing in the limit e — > 0. With a little manipulation of the summations becomes, 

P» = 2 (f + loge)+/ 1 - £ da;x z - 1 + ^ + T ^:-2E: = o^ u 9z a 

It is worth pointing out that in dropping the terms vanishing in the limit e — > the cut off e 
has now disappeared altogether from the calculation. The Euler function is defined as 

il)(x) = — logr(x) (4.25) 

dx 



( 7 )The following dilogarithm identity is required, Li 2 (x) = ^ — log (x) log (1 — x) — Li 2 (1 — x) to obtain the 
term in the form of jSj. 
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hence, 

V> (x + 1) - if> (x) = - (4.26) 

x 

(r (x) is the Gamma function). In terms of the Euler function, 

p» = (z + 2) - i) (z) + V (3) + V (1) ■ (4-27) 
Using the identity 14.261 again and expanding T (x) around x = 1 we find 

p» = -ij) ( z + 2) - 4j (z) + - - 2 7 . (4.28) 

Finally, substituting 14.281 into 14.22] we have an explicit expression for the structure function in 
terms of its Mellin transform (|4.19jl 

D NS (x,rj) = — exp(^-- 7 JJ dz x~ z exp [-7, (z + 2) - $ (z)) J . (4.29) 

An analytic expression for the integral above is not known. The integral may be performed in 
the soft limit x — > 1, this is known as the Gribov approximation. The factor a; -2 shows that 
the integral is dominated by large values of z. The Euler Gamma function is defined, 

T(z) = (z-1)\ (4.30) 

in the large z limit one can approximate the factorial operation by Stirling's formula n\ pa 
\^2tvz z+ 2 exp (—z). In the large z limit we have ip (z) ~ logz 

^^(,,,)=-exp(-,( i - 7 ))/_ d2 OT. (4.31) 
By substituting y = zlogx, 

1 /1/3 \\ 1 1 rc'—ioo j 

(x, r/) = — exp -r/ - - 7 (logx)^" 1 / dy ^ exp (-y) . (4.32) 

After the transformation the limits have changed c + ioo — > clogx + ioologx = c' — zoo. 
Importantly c 7 is negative as logx < and c > 2 > is required for the Mellin transform to 
exist. The integrand has two obvious singularities, one at y = and another at Re (y) — > —00 
(Re (2) — > 00). For \y\ — > 00 and Re(y) 7^ —00 the integrand is zero, this being the case we 
can close the integration contour by joining up the two ends at d + ioo and d — ioo such that 
it becomes a hemisphere of radius \y\ = 00 in the Re (y) > d region of the complex y— plane, 
as the integrand is zero along this addition to the contour. Given the contour is closed in the 
complex plane we can deform it as we please so long has we keep all poles inside it as its value 
is given by lui times the residue at the pole, in this case the only pole is at y = 0. The Hankel 
contour C-^ in the complex plane is an open contour which comes in along just under the real 
axis from +00 goes around the origin and back out to —00 just above the real axis. We can 
clearly deform our contour to Hankel's contour as it contains the origin and no other poles and 
because our integrand is zero at \y\ = 00 and Re (y) 7^ —00 we can open the contour again at 
\y\ = 00 making it exactly the Hankel contour. The integral can then be brought into the form 
of Hankel's integral representation of the Gamma function 

r&y^—Jiyi-yY^^-y)- (4.33) 
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Inserting 1 = (_i)+i r '- 1 (-1) +1 and EM into EM gives 



D^ lbov (s, rj) = exp ( -77 ( - A - j ) ) v ~ 7 N • (4.34) 



2 ' V4 



(1-x 



(we have used lim^! — log (x) = linx^! 1 — x). 

The result of the Gribov approximation constitutes the perturbative result derived earlier 
extended to all orders with the caveat that it is only valid for the limit x — > 1. To quote 
Gribov and Lipatov [IH] the approximation above "can be considered as a generalization of 
the Sudakov form factor". The perturbative structure function calculated earlier has no such 
constraint on it but is nonetheless a finite order perturbative result. Ideally one wants a single 
expression for the structure function which tends to the Gribov result at in the soft limit and 
to the perturbative result away from x —> 1, the desired expression will interpolate between the 
perturbative structure function and the Gribov approximation. Such interpolation constitutes 
what is known as exponentiation. Integrating 14.341 over the soft phase space gives, 



dx D^ ibov (x,rj) = exp ( - V [ - - 7) J ^ . (4.35) 



We denote this Dg^ t (e, 77) and rewrite the e term exp ( | log e 



exp I -±777 I ( \ /3 
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Dsoft V) = / VV exp ( ^ ( ^ + loge ) ) . (4.36) 

r [l + - 2 t) 

The soft part of the structure function Dg^ t (e, 77) gives an expansion in powers of r\ with 
coefficients exactly the same as the coefficients of 77 in the delta term obtained earlier by 
iterating the Lipatov equation, so justifying the earlier quote of Gribov and Lipatov. We have 
checked this explicitly to 0(/3 3 ). Had we calculated the O (/5 4 ) delta term with Lipatov's 
equation we would see its coefficient is the same as the coefficient of the O ((3 A ) term in the 
expansion of Dg^ t (e, 77). In the perturbative solution the delta terms give the structure function 
in the 1 — e < x < 1 region. If we were to integrate the delta terms over 1 — e < x < 1 and 
expand in 77 we would get the same result as we would just expanding Dg^ t (e, 77) up to the 
same power in 77. Denoting the delta term in the perturbative expansion calculated to O (rj 1 ) 
as 6 (1 — x) Aw (e, 77) we have, 

f dx 5 (1 - x) A« (6, 77) - f dx D^ bov (x, 77) ~ O ( V t+1 ) (4.37) 

this suggests that we modify the structure function by hand to include higher order soft effects 
by replacing the delta term altogether with Dcribov (, x > v) > 

5(l-x)A^(e, V )^D^ bov (x,r 1 ). (4.38) 

How does this affect the theta term etcl With this replacement we can safely take the already 
implied limit e — > 0, the theta function for the theta term is then replaced by 1. Then the most 
obvious thing to do is to demand that, up to the same order in 77, the new structure function 
is the same as the old purely perturbative one was in the region < x < 1 — e and not care 
about the terms that are higher order in 77. Expanding in 77 we have 



^Gribov ( X iV) — 1-x (4) 



3+41og(l-j!) 



(r) 

1-x V 4 / 

+ 

+ 0(77 4 ) 



f-|7r 2 +61og(l-x)+41og 2 (l-x) 



(4.39) 
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Recall that the Gribov solution is valid in the large x limit. The perturbative solution obtained 
order by order in r\ by iterating solutions through Lipatov's equation requires no such approxi- 
mation. Consequently one expects that the perturbative solution and the Gribov approximation 
should agree in the limit x — > 1. This is indeed the case, it is easy to verify that the x — > 1 
limit of the theta terms obtained in the perturbative case are equal to those obtained in the 
non-perturbative case. If we denote the ith order in rj of a quantity by appending superscript 
(i) to it we can define an improved structure function 



^(0 (a . v) = { D oJiL (*, V) + (d NS{ * ] (*, V)\ e _ term - D*™ (x, r,)) < i < 3 



(4.40) 



i.e. 



D' NS (x, V ) = D™ bov (x, V ) + E (*, »/) | term - dS£ ^)) • (4-4i) 

i=0 v J 

This is the exponentiation prescription of Kuraev and Fadin The condition on i relates to 
the order to which the perturbative solution was found, in our case 3. Applying prescription 
14.401 to the perturbative theta terms obtained so far is trivial albeit tedious subtraction, we 
find agreement with the results in [EllZj. It is of note that the integral of this function over the 
range 1 — e — >• e is finite in the limit e — > 0. The integral is finite because, as stated above, in 
the divergent x — > 1 limit the theta terms are equal to the corresponding Gribov terms inside 
the summation, in addition the integral over the all orders Gribov approximation to the left of 
the sum is finite (see equation I4.36|) . 

An alternative exponentiation prescription is that of Jadach and Ward called YFS exponen- 
tiation. Here the only difference is that the Gribov term is extracted from the iterative solution 
as a factor rather than subtracted as in the Kuraev Fadin exponentiation viz 

n'NS(i) (rr „\ _ / £j+fc=i D Gribov (x, V) A (fc) (x, rj) < i < 3 f 

^'iD^M^Hx,,) i>3 ' (442) 

where A^ fc ^ (x, rj) is defined by 



D NS{i) (x, V )\ e _ term = E D^(x,r))AW(x,r)) (1 < j + fc < 3, < k < 2) . (4.43) 

j+k=i 

This is a system of linear equations that is easily solved for the A( fe ) (x, 77), 
A(°)(x,77) = 

am(x, V ) = "";;x7!' 8 r m - m (4.44) 

D„ .) Ix.ti) D„ .; (x,ri) 

1 ' " ~~ D^^fxr,) Z3^ S <S 'far?) [ 1 1 ) D NSW (x V ) 1 

Gribov^ ' Gribov^ ' Gribov^- ' '-' 

giving 

D /7VS (x, 77) 

= (x, 77) x (§ (1 + x) 2 + i (-2 (1 - x) 2 - (1 + 3x 2 ) logx) 77 

+ 3^ (l2 (1 - x) 2 + (6 (1 - 4x + 3x 2 ) + (1 + 7x 2 ) logx) logx - 12 (1 - x) 2 Li 2 (1 - x)) rj 2 ' 

(4.45) 

This is the result given in jHllZj. 

To summarize, we have re-derived in this section again the structure functions of Kuraev 
and Fadin l4~4"0l and Jadach and Skrzypek. At no point is the underlying hard scattering process 



D ris 0){x,rf} 


6 — term 


D" S ^ ( 

Gribov ^ 
D NS ^{X,7 1 ) 


e-term U Gribov*^ > 


D™.W ( 

Gribov * 

D NS ^(x, v ) 


D GriboJ X ^) 
8-term ^ Gribov^ >'l > 


D NS(1) (x n) D NS(1) (x n) 

Gribov^- ' Gribov^- ' 
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Radiator Function 



Figure 2: The plot above shows the radiator functions R\ (s, s') with s = (200 GeV) 2 , obtained 
from the ZFITTER package by extracting it separately for each angular bin just as in figure d 
but this time including initial state-final state interference effects. 

referred to, with this approach folding the radiator function with New Physics is just as valid 
as folding it with the standard model. In addition, within the context of Lipatov's equation 
(i.e. for the case of near collinear emission of radiation) the calculations suggest that the 
factorization of the corrections holds at the level of the differential cross section not just the 
total cross section. We are seeing this effect in the flatness of the radiator functions we have 
derived in figure [TJ In fact the third order Kuraev-Fadin structure function 14.401 is folded with 
itself 13.51 and present within ZFITTER as the default setting for the flag higher F0T2 which 
governs the implementation of higher order radiative corrections. Crucially at this level of 
approximation the hard scattering interaction is not referred to, this is clearly good news for 
fits to new physics. Approximations are made in the leading log approximations e.g. the finite 
order perturbative result nevertheless they have been found to agree very well (KF better than 
0.2%, YFS better than ~0.01%) with exact numerical solutions of the Lipatov equation. This 
would seem to vindicate the use of ISR corrections, the numerical radiator function, derived 
from ZFITTER and other such packages. 

5 The Breakdown of the Simple LLA Structure Function 
Approach. 

In general the standard model predictions obtained from ZFITTER and the other standard model 
difermion packages allow many more sophisticated radiative corrections to be applied to the 
process in question besides just ISR. In particular initial state-final state interference (ISR*FSR) 
effects corresponding to interference between diagrams with bremsstrahlung emitted from initial 
state particles and diagrams with bremsstrahlung emitted from final state particles, as well as 
box diagrams resulting from emission and re-absorption of photons, have a significant (~10%) 
effect at the level of the differential cross section despite having only a negligible (<1%) effect 
on the total cross section. The radiator functions derived from ZFITTER with these corrections 
included are shown in figure [2] where they have acquired a significant angular dependence that 
was not seen in figure [0 

This angular dependence would seem to ruin the argument put forward in [J which says that 
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these corrections should be independent of the polar angle if in fact they are independent of the 
hard physics process. It also highlights the weak point in the leading log approximation, the 
leading log approximation starts to break down for exclusive observables i. e. highly differential 
quantities and measurements obtained with strong cuts pj. ISR*FSR interference corrections 
clearly fall outside the mandate of the structure function / leading log approach which describes 
ISR effects so well. 

In fact the cross section obtained with angular cuts requires two radiator functions. The cos 6 
symmetric and cos 6 asymmetric parts of the tree level cross section are convoluted separately 
with different radiator functions and the result added [TH1 EH EH- The cross section in an 
angular bin i with edges at c = cos6> = q, q+i o~i (s, Cj, Ci+i) is 



Ci+1 d < ^°" 
i dc' 



Defining the symmetric and asymmetric cross sections cr^ (s, 
tively in the region — q < c< q as 



(5.1) 

-Cj, Cj) and &fb (s, — Cj, q) respec- 



a T (s, -Ci,Ci) 
o-fb (s, -Q,q) 

the cross section in the angular bin is 



Jo* ^ + J_° Ci dc J 
Jo" dc f - J° 



dcf 



(5.2) 



ctt (s, -Ci+i, Q+i) + (Tfb (s, -Ci+i, Ci+i) - cr T (s, -q, - CTfb (s, -Cj, Cj) . (5.3) 



In terms of a cross section corrected for radiative effects we have 



a T [s, -Ci, Ci) 
a FB (s, -Ci, c^ 



/// ds' R T (s, S', Ci) CTree (s') 

min 

//, . ds' .Rfb (S, S', Ci) (J T ree (s') 



So far we have assumed 

i? T (s, s', Ci) = R FB (s, s', c^ , 
and that to good approximation (|3.8j) . 



(T (S, Cj, Cj+i y 



ds' i?i (S, S') (TTree (s', Cj, C i+1 ] 



(5.4) 
(5.5) 
(5.6) 



where, recalling section [31 i? (s, s') is a universally applicable radiator function. According to 
Bardin et al [11] "Strictly speaking an ansatz like 15.51 is wrong." In [TJ] the authors go on to 
say that near the Z peak the ansatz 15.51 nonetheless gives "excellent agreement with the correct 
result". The authors say that this effect is due to the dominant corrections being from soft 
photons. They go on to say that if an s' cut is in use the agreement is even better. We appeal 
to this soft photon dominance ansatz due to the restrictive s' cuts used in obtaining samples 
to fit with, \fs' > 0.85a/s, the degree to which it holds is born out by the flatness in cos^ of 
the corrections we derived, R (s, s') in figured If the corrections were significantly different for 
the cos 6 symmetric and asymmetric parts of the differential cross section R (s, s') would not be 
flat in cos 6 as we essentially are just dividing the corrected cross section in each angular bin 
by the born cross section in each angular bin, bin by bin in s', the born cross section could not 
possibly factor out. Prior to considering initial state-final state interference effects there seems 
to be nothing wrong in using the same radiator function for c even and c odd parts of the cross 
section, the structure functions we derived describe the radiator function very well. 

As you can see in figure El the addition of initial state-final state interference effects seem 
to ruin all our previous hypotheses. Ideally one would hope that just because the radiative 
corrections have acquired this angular dependence they are still nonetheless essentially universal 
like the other (ISR) corrections or at least largely process independent in our phase space. Given 
that we are dealing with box diagrams such factorization seems like wishful thinking. 
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Figure 3: ISR*FSR terms in the cross section originate from diagrams with photons emitted 
from the initial legs and final state legs interfering. 




e" f e" f 



Figure 4: In addition to the bremsstrahlung diagrams in figure 01 at the same order in a, 
ISR*FSR terms in the cross section also originate from diagrams with photons box diagrams 
interfering with the lowest order diagram. 

6 ISR*FSR Universality. 

In this section we describe how the ISR*FSR interference corrections, despite giving the naive 
radiator function an important angular dependence, is nonetheless mostly a universal correc- 
tion. The ISR*FSR interference corrections for some underlying hard scattering involving the 
exchange of an unknown particle (???) may be represented at O (a 3 ) by the Feynman diagrams 
in figures 01 and HJ 

In subsections 16.11 and 16 . 21 we will use the Feynman rules of [2] and the kinematic notations 
of appendix A. In these subsections we will try to maintain an explicit notation. The matrix 
elements of the two fermion events under discussion at this order have four fermions and 
hence four spinors in them, v e (p + ) , u e (p_) represent the incoming e + and e~ , Vf (q + ) , Uf 
represent the outgoing antifermion and fermion respectively. Given that we are trying to discuss 
radiative corrections in such a way as to find and highlight any universal features of them it 
will therefore be useful to denote the born level process generically as 

A B om = v e (p + ) ...Hard...i> f (q + ) (6.1) 

where ...Hard... is defined to contain the electron and fermion spinors and more importantly, the 
propagator of the exchanged particle and its couplings to the incoming and outgoing particles 
(see also equation 16 .3j) . For example in the case of QED we have simply 

...Hard... = iQ e e^u e (p_) % ® >w 2 u f (g_) iQ f eY- (6.2) 

(P- +P+) 

6.1 O (a 3 ) Soft Bremsstrahlung. 

The contribution of soft bremsstrahlung corrections to 2-fermion processes are greatly impor- 
tant. By soft bremsstrahlung we mean bremsstrahlung which carries off an amount of energy 



15 




Figure 5: Some O (a 3 ) bremsstrahlung correction diagrams. 

smaller than the energy resolution of the detector, soft photons carry off a negligible amount 
of energy. Given that the detector cannot therefore tell the difference between a genuine O [a) 
process and an O (a 3 ) soft bremsstrahlung correction to it, we must add such contributions to 
our predicted cross section. The amplitude for a general difermion process with bremsstrahlung 
emitted from the external positron leg of the diagram (figure EJ) will be of the form 

A Brem = e* x (k) v e (p + ) {zQ e e 7 A ' ( 7 (i6+ ~ f ) + 7 e) | ...Hard..M f (q + ) , (6.3) 

with Q e = — 1 for an electron (generally Qf is the charge on the fermion/antifermion / in units 
of e). If we now work in the limit that \k\ — > 0, s' — > s, the soft photon approximation, we see 
that 

Asoft = -Q e ee* x (k) v e (p+) 7 A ^r ...Hard...v f (q + ) . (6.4) 

Simple 7 matrix manipulation and the application of the Dirac equation to 16.41 gives 

(6.5) 



A*oft = -2tk^Qe ee \ ( k ) v * (P+) ( 2 P+) -Hard...u f (q+) 

Qe& f~ p + Ab orn 



From this we see that in the soft photon approximation the bremsstrahlung corrections to 
the born level process are independent of what that born process is, we have extracted the 
bremsstrahlung correction as an over all constant factor while staying blissfully ignorant of 
the underlying process. Hopefully it is clear that these tricks can be employed with the other 
external legs in considering bremsstrahlung emitted from them in the soft photon approximation 
giving for bremsstrahlung emitted from the ith external leg of the diagram 



A^Soft — Qi e 7 Abotu (6.6) 
K.pi 



where pi denotes the momentum of the external leg. Consequently, taking into account the 
bremsstrahlung emitted from all different legs we find the matrix element 

AsoftAsoft = QiQje 2 y^*^ ^ ^ A ] Born A B om- (6.7) 

Summing fermion spins and polarizations of the photon (j2 po is e^e* u = —g^ and inserting the 
phase space and flux factors factors gives in the k m soft photon approximation 

ct d 3 fc I 

da = ~W § QlQm ' P ^ \k. Pl ) (k.p/ aBorn - (6 - 8) 

The soft photon cross section is the contribution to the cross section from soft photons where soft 
photons were defined earlier as being those low enough in energy to render the bremsstrahlung 
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process indistinguishable from the tree level process, we denote this energy u, the soft photon 
cut-off . Integrating over the final state phase space of the soft photon we have 



a (if d 3 k 1 \ 
dasoft = Y,QiQjPi-Pj \T~ /i-i T~Ti — \Ti V da Bo m- (6.9) 

Note that this differential cross section contains contributions from pairs of diagrams which 
constitute initial state-final state interference. This would appear to be the desired result, 
the radiative corrections are factorized from the hard scattering process at the level of the 
differential cross section. 

On its own the result 16.91 arguably does not make much sense as the phase space integral 

1 f d 3 k 1 , 

(6.10) 



Atx J\k\<u> k (k.pi)(k.pj) 

is infrared divergent, as k —> dcr — > oo . The infrared divergence of the ISR*FSR bremsstrahlung 
contribution actually cancels an infrared divergence from the ISR*FSR box diagrams. The 
phase space integral must be expressed in a form whereby we can add the bremsstrahlung and 
box contributions to cancel the divergences before we can fully assess the universality of the 
ISR*FSR interference. We have treated the divergence with dimensional regularization. We 
replace all four dimensional dependence with the usual prescription 4 — > n = A + e' 

1 ' ^ 1 W * , (6.11) 



An J\k\<u> k (k.pi) (k.pj) An J\k\<u> k (k.pi) (k.pj) 

The factor p, N is introduced to keep the overall mass dimensions of the integral the same as 

N i.e. N= -e' 

(6.12) 



they were before extending to N dimensions. p N has mass dimensions 
i.e. 16.111 is 

1 flY' r d e ' +3 k 1 



An \n ) J\k\<w k (k.pi) (k.pj) 

This photon phase space integral is technical and not so illuminating for our present discussion 
so we shall confine it to appendix B. After completing the integral the contribution to the 
differential cross section from O (a 3 ) soft bremsstrahlung diagrams is found to be 



2a 

do~SoftJnt — 

7T 



da Born- (6.13) 



This is the universal soft photon bremsstrahlung correction to the differential cross section 
from ISR*FSR interference. Left as it is above, the correction to the differential cross section 
is infrared divergent, there are terms O (e /_1 ) and the limit e' — > must be taken to return to 
four dimensional Minkowski space. For our purposes the essential feature of this result is that 
the radiative correction is factorized from / independent of the new physics exchange process. 



6.2 O (a 3 ) Virtual Diagrams. 

The infrared divergences (jj in the soft bremsstrahlung result are well known to be physical, 
that is to say that they are a real part of the theory as opposed to the ultraviolet divergences en- 
countered in renormalization which are a consequence of our ignorance of short distance physics, 
likewise the infrared divergences cannot be subtracted ad hoc. The infrared divergences due 
to the emission of soft photons do however cancel against the corresponding virtual diagrams, 
it is a well known fact that the soft photon contribution to the cross section from initial state 
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Figure 6: The box diagrams contributing to the ISR*FSR interference cross section. On the left 
hand side we have direct box diagrams, from top to bottom these are represented by amplitudes 
Adi and Am respectively. On the right hand side we have crossed box diagrams, from top to 
bottom these are represented by amplitudes Aci and Ac2 respectively. 

bremsstrahlung cancels against a corresponding infrared divergence from the corresponding 
vertex diagram which involves emitting a photon from one initial state leg and absorbing it 
on the other. Likewise the initial state-final state interference cross section arising from the 
interference of initial and final state bremsstrahlung diagrams does not make much sense on 
its own, it must be added to the corresponding virtual graph. The virtual graph for ISR*FSR 
bremsstrahlung is by analogy to the case of initial state bremsstrahlung the one formed by 
taking the photon emitted from an initial state fermion and joining it to a final state fermion. 
This gives the two diagrams of figure H3 

At O (a 3 ) these diagrams only contribute by their interference with the tree level (hard) process. 
Unlike the soft bremsstrahlung corrections the ISR*FSR box diagrams do not simply represent 
multiplicative corrections to the born level amplitude. These diagrams were studied in the 
the hope that, though they do not fall under the banner of universal radiative corrections, 
they may have some component which is universal. It is quite easy to show that the box 
diagrams have universal, factorizable corrections. The tricks previously applied to the soft 



18 



bremsstrahlung diagrams also work, to some degree, with box diagrams (provided there is an 
internal photon). Let us write down the amplitude for the first direct box diagram in figure 
were we feign ignorance of the couplings and propagator of the X° boson^ 8 ). We denote 
the couplings associated of X° some particle A as Ca and propagator contributions to the 
numerator as Vx°- We assume the denominator of the propagator is the usual {k 2 — m 2 x ) form 
(for an X° particle four momentum k and mass nix)- The couplings can be considered to 
have Dirac and Lorentz indices, the Dirac indices are contracted with the rest of the Dirac 
algebra in the numerator and the Lorentz indices are contracted with those in Vxo. With these 
notations Vxo---C e u (p_) u (p + ) Cf.../ (s 2 — m 2 x ) corresponds to ...Hard... in the case of the soft 
bremsstrahlung. Finally we shall assume the interaction is charge conjugation invariant i.e. 
Ca = —Ca- The amplitudes for the diagrams in figure El are then 



Am = 


-ie 2 j 


d 4 k 

(2^) 4 




A D2 = 


-*e 2 f 


d 4 fc 
(2^) 4 


V x o 


Acx = 


+ie 2 f 


d 4 fc 
(2n) 4 


V X o 


A C 2 = 


+ie 2 f 


d 4 k 
(2n) A 


V X o 



v(p+h^(-¥-rf++™e)C e -u{p^)u(q^)Cf( y -k'-tf + +m f )-y ll v(q + ) 



k 2 (k 2 +2k.p + ){k 2 +2k.q + )(k 2 +s+2k.(p++p-)- 



(6.14) 



u(p+)C e+ (^-+m e )^u(p-)u{q-) 1 ^(^-+m f )C f u(q + ) 

k 2 {k 2 +2k.p-)(k 2 +2k.q-)(k 2 +s+2k.{p + +p-)-m 2 x ) 
9(p+)C e+ (k'+p'-+ me ) 1 »u(p-)u(q_)C f (-y-q , + +m f ) lll v(q + ) ' 

k 2 (k 2 +2k.p-)(k 2 +2k.q + )(k 2 +s+2k.{p + +p-)-m 2 x ) 
9(p+h»{-k'-tf + + me )C e -u(p )u(q-) 1 ^(k'+4-+m f )C f -v{q + ) 
k 2 {k 2 +2k.p + )(k 2 +2k.q-)(k 2 +s+2k\p + +p-)-m 2 x ) 

where we have defined the loop momentum in each case such that it flows counter-clockwise 
with the X° propagator carrying momentum k + p + +p-i. e. the propagator is the same in each 
diagram. We have defined the loop momentum k such that in each diagram it is the photon 
four momentum flowing into the incoming antifermion line. Firstly condense the numerator as 
much as possible. The terms in the numerator may be quickly got rid of by anticommuting 
them so that the on-shell Dirac equation may be used. Considering the left hand side of the 
numerator of Adi 

v{p+)l^{-fc- i>+ + m e )... ( 
= eU{p + ){-Y/k + 2^ + )... ■ [0 - L °> 

Using the same techniques on the right hand side of the various numerators the amplitude can 
be simplified to 

HP+)(-y"¥+2p^)C e -u{ P -)u('}-)Cf(¥^+2q+„)u(q + ) 
k 2 {k 2 +2k.p + )(k 2 +2k.q + ){k 2 +s+2k.{p + +p^)-m 2 x ) 
P(p+)C e -(¥^+2p l i)u(p_)u(q_)('y^+2q.„)C f u(q + ) 
k 2 {k 2 +2k.p_)(k 2 +2k.q-)(k 2 +s+2k.(ji + +p-)-m 2 x ) 
9(j> + )C e -(^+2p»_)u(j,-)u(q_)C f (¥~ /fl +2q +fi )v(q + ) ' 
k 2 (k 2 +2k.p-)(k 2 +2k.q + )(k 2 +s+2k.(p + +p-)-m 2 x ) 

k 2 (k 2 +2k.p + ){k 2 +2k.q^)(k 2 +s+2k.{p + +p.)-m 2 x ) 

Assuming that the couplings (C e //) and numerator of the X° propagator Vxo do not depend 
on the loop momentum we can decompose the amplitudes into tensor (T Q/ g), vector (V a ) and 
scalar (S) loop integrals according to the number of powers of k in the numerator 
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= -*e 2 f 
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(2nf 
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= +*e 2 J 


d 4 k 
(2^) 4 


V X o 


A C2 


= +*e 2 f 


d 4 fc 
(2-) 4 


V x o 



(6.16) 



A D1 = -ie 2 u(p + )^ a ...-f l3 - ff ,u(q + ) xT aP 

-ie 2 v (p+) ( 7 ^...2g +M + 2p%...^ lp ) v (q + ) x V a , (6.17) 
-Aie 2 (p+.q+) v (p+) ...v (q+) x S 



(8)p or now i e t us assume that X° has the generic ™ high energy boson behavior i.e. the diagram is UV finite. 
We discuss this point more later. 
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where 

f d 4 fc 

V <* = J (a*) 4 fc 2 f f fc+P4-) 2 -m2 1 f (fc+o+) 2 -m 2 1 f (k+v++V-? -ml ) ( 6 - 18 ) 
S = J 

and C e -w (p_) u C/ has been replaced with for brevity. We can see that the scalar 
term in the amplitude of 16.171 is a factor multiplied by the born scattering amplitude of the 
X°particle 

2it x e 2 x S x (s - m x ) x ^P + ) C e -u (p^ u( q ^) C f v ( q+ ) ^ = 2itx # xSx ( a _ ^ x ^ 

(6.19) 

The numerator algebra has given a universal^ 9 ) factor which multiplies the tree level numerator 
algebra 4(p + .q + ) = —2t however the same is not true of the denominator. The factor S x 
(s — m x ) depends on the underlying hard scattering process through the dependence of S on 
m x , the mass of the exchanged particle. Using the on mass shell relations we can rewrite the 
denominator of the integrals in the |A;| — > limit as 



llm |fc|-0 fc2 ^ fc+p+) 2_ mi ^ (fc+9+) 2_ m 2)( (fe+p++p _ ) 2_ m ^) ^ 



k 2 {2k.p + ){2k.q + )( y s-m 2 x ) 

By trivially counting powers of \k\ we see that in limit \k\ ^ the tensor integral goes as 
/ d \ k\ \k\, the vector integral goes as / d \ k\ and the scalar integral as / d\k\ \k\~ . Consequently 
the scalar integral is infrared divergent and the other integrals are infrared finite. Simple 
power counting also shows that all of the integrals are finite in the ultraviolet limit \k\ — > oo. 
Expanding the last bracket in the denominator we can write, 

r r\ 4 k 1 

S= / ————j ~, r-, r . (6.21) 

J (2vr) 4 k 2 ((Jfe + p+f - ml) [(k + q+f - mfj (k 2 + s + 2k. (p+ + p_) - m 2 x ) 
In the infrared divergent \k\ — > limit we can rewrite the last part of the denominator 

V + s + 2k. (p+ + p_) - m 2 x ) = (s- m\) , (6.22) 
therefore the radiative correction factor in the scalar integral term 

lit x e 2 x S x (s - m x ) , (6.23) 



has a divergent piece which is universal because the constant factor (s — m 2 x ) in 16.231 cancels 
that which appears in the divergent limit of !6.211 We can perhaps express this result in a better 
way writing 



m 



x, 



= 1 fc 2 + 2fc. (p ++P .) 

(k 2 + 2k. (p+ + p-) + s - m 2 x ) {k 2 + 2k.(jp+ + pJ)+s-m 2 x y V ' ; 

which makes the scalar integral term equal to 

2itxe 2 x S x(s- m\) = J ( ^ 4 fe2 ( (fc+p+) 2 J^ (fe+9+)2 _ m 2) 

f d 4 fc 2itxe 2 x(k 2 +2k.{p++p-)) ■ (6.25) 

~~ * (2tt) 4 k 2 ({k+p + ) 2 -m 2 )((k+q + ) 2 - m 2 )(k 2 +s+2k.(jp + +p„)-m 2 x ) 



(9 ) Universal in the sense that we have not referred to the specific details of the X° exchange process. The 
amplitudes corresponding to the other diagrams (Ad2 , Ac\ , Ac2) also give such factors the only difference 
being the replacement t <-> u for crossed box diagrams i. e. the differences in the universal factors are due to the 
topology of the diagrams rather than their physics. 



20 



The first term on the right hand side of 16 .251 is universal, the terms in the denominator corre- 
spond to the photon propagator and two fermion propagators and the —It in the numerator 
came from the universal part of the numerator algebra. The second term on the right hand 
side of 16.251 clearly depends on the hard scattering process (note the m x terms), it is actually 
a vector and a tensor loop integral. The universal term goes as / d \k\ \k\~ l in the \k\ — > limit 
while the other term goes as / d \ k\ i.e. only the universal part of the scalar integral is infrared 
divergent. The degree to which the box diagram corrections are universal depends on how big 
the first integral is in 16.251 relative to the other terms. If we denote the universal part of the 
scalar integral coming from the denominator of the amplitude 

Su q+) = I ^k^((k +P+ f- m A((k + q+ f-^y (6 - 26) 



m 



we can then rewrite the amplitude 16.171 



Am = -ie 2 v (p+) YT-1 P 1^ (<?+) T aP - 2ie 2 tT«A Bo m 

-ie 2 u (p+) (j^ a ...2q +fl + 2^... 7 a 7/ J v (q+) V a - AieHV. (p + + p_) A Born . (6.27) 
+2ieHS u (p + ,q + )A B 

orn 

where represents C e -u{p_)u(q_)Cf. Hopefully it is clear that the decomposition of such 
scalar integrals into a universal part and non-universal part is as general as propagators with 
denominators of the form p 2 — m 2 , that is to say no matter what X°is we will always find a 
term Su 16.261 

We can repeat this process with the other diagrams (figure EJ) in exactly the same way, in 
each case we have the generic result that the radiative correction from ISR*FSR box diagrams 
factorizes from the tree level process in the infrared divergent part of the scalar loop integral. 
Referring back to 16.161 we see that the other topologies of the process give the following scalar 
terms 

-A-Dllscalar = +^ 2 tS V (p+, q+) A B orn 

A D2\ Scalar = +2ie 2 tSu{p_,q.)A Born < 628 

■Acilscalar = -^uS V (p_, q+) A Bo rn ' 
• A C2\ S calar = -^uSu (p + , g_) A Born 

Again we have found what we are looking for, factorization of the radiative corrections. The 
ideal result, namely that the differential cross section due to ISR*FSR box diagrams interfering 
with the tree level process is of the form 

dcr Bo x = Universal Factor x da Born , (6.29) 

i.e. exactly factorizable requires that the vector and tensor integral terms as well as the non- 
universal parts of S are negligible relative to the infrared divergent, universal piece of S. Clearly 
an exact factorization is impossible, factorization will only ever occur to some degree. The 
question then becomes, is a good degree of factorization possible and generic? Requiring a 
universally good degree of factorization means that the factorizable part of the amplitude must 
somehow dominate all of the other parts. It does not seem improbable that this be the case as 
the universal correction we are interested in always corresponds to an infrared divergence, in 
fact it corresponds to the only divergence, infrared or otherwise. We shall briefly postpone the 
discussion regarding the relative sizes of the various contributions to da Box to discuss an initial 
simplifying assumption. 

It is important to note that some of the new physics X° particles which we hope to apply 
this analysis to will have momentum dependent couplings e.g. gravitons (gravitons couple to 
the energy momentum tensor). Another, more pertinent point is that bosonic propagators 
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generically contribute polynomials in the momentum of the particle they represent to the nu- 
merator as well as the denominator of the amplitude. This does not change our current analysis 
except at the point of decomposition into scalar vector and tensor integrals 16.171 as gravitons 
will give C e -, Cf and Vx° then have a dependence on the loop momentum. In this case the 
same decomposition can be done into scalar, vector and tensor integrals and crucially one finds 
that the scalar term is the same as in equation 16.271 To see this all one has to do is simply 
multiply out any k dependence in Vxo , C e -, Cf i.e. decompose Vx°, C e -, Cf into scalar, vector 
bits etc. Another way to understand this is to recall that the scalar part is infrared divergent 
and the momentum in the photon in our diagrams is just k the loop momentum, so naively 
one can think of the photon line in the diagrams as vanishing in the infrared divergent limit. 
Thus even in the case where there is a non-trivial dependence of the couplings and propagator 
of X° on the loop momentum the generic infrared term still occurs, essentially one just has to 
expand in the loop momentum. In the case of momentum dependence of the couplings etc the 
tensor integrals will naturally involve tensors of higher rank than just two i.e. we will have 
tensor integrals of the form 

d k kgkpky... , , 

(2ix) k 2 \{k + p + ) -m 2 .)[(k + q + ) - mjj ((k + p + + pS) - m x J 

This causes us to rethink the divergent structure of the amplitude, up to now only the term 
involving the scalar integral, the term that represents the universal radiative correction was 
(infrared) divergent, all other terms were ultraviolet and infrared finite. This was encouraging 
from the point of view that this could be expected to make the universal correction the dominant 
one. In the most general scenario momentum dependent couplings and the numerator of the 
propagator could give ultraviolet terms. Power counting shows that one requires four powers 
of the loop momentum in the numerator to have an ultraviolet divergence. In Giudice et al 
[28] the Feynman rules are derived in the unitary gauge in which the graviton has a propagator 
similar to that of the standard model gauge bosons in unitary gauge, that is to say the high 
energy behavior goes as . We shall suppose that in this model (or at least in a fully consistent 
theory of gravity) it is possible to work in a gauge analogous to the the gauges in which 
the graviton propagator goes as \ at high energy. With this assumption then we have at most 
four powers of the loop momentum in the numerator of any loop integral, two from the internal 
fermions and one from each coupling of the graviton to the fermions. Such integrals diverge 
logarithmically. Assuming we cut the integral off at \k\ = A they should contribute logarithms 
~ log — (A and ^fs are essentially the only two energy scales in the diagram) . The model of 
|2E| has a cut off of order 1 TeV. For A = 10 TeV, s = 200 GeV we have log ^ = 5, this is 
something we should be mindful of when considering the size of any infrared divergences. 

The infrared divergent universal term must cancel the corresponding infrared divergence 
from the bremsstrahlung calculated in the last section. To see this we have to regularize the 
integral. We extend the dimensionality of the integral from 4 to n = 4 + e' dimensions as in 
section 16.11 

Su (p + , q + ) = -i-j lim (-) [ d i+ *'k —j- -y. (6.31) 



k 2 ((k + p + ) 2 - mfj ((& + q + ) 2 - mj^j 



The denominators can be combined through the usual method of introducing Feynman param- 
eters 



1 f 1 f 1 

k 2 [{k + p + f - ml) ({k + q + f - m 2 f ) Jo Jo 



1 



o " (k 2 + 2yk. (xp+ + (1 - x) q + )f 

(6.32) 
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To use the standard integrals we complete the square in the denominator and shift the variable 
of integration by a constant amount, fc M — > k 1 * = + y (xp+ + (1 — x) g+) to give 



S v (p+,q+) 



-in 2 , I 1 
hm 



(2vr) 4 e'^o \n) Jo 



1 dx f 1 dy y tt^'T ( 1 - \e' 



1 



^y 2 (xp+ + (1 - x) q+Y 



(6.33) 



From now on we denote A (p + , q + ) = (xp + + (1 — x) q+) 2 . In the current form we can perform 
the integral over y, 



1 2 



€ 



(6.34) 



Substituting this into Sjj (p+, q + ) and expanding the term in brackets about e' = gives 



Su q+) = ~ 32^ &o f dx 1 (? " bg (t) + bg (f 
Expanding A (p + , q + ) (keeping all masses) 

A (p + , q + ) = (x 2 t - x (t - m 2 + m 2 A + mfj , 



(6.35) 



(6.36) 



we see that Su (p+, q+) is really a function of the Mandelstam variable t. Likewise A (p + , q + ) 
is really a function of t, A (t). We can perform the integral by completing the square in A (t) 



where 



A(p+,q + ) = -t (k 2 -y 2 ) 
y = x -|(l + 5) 
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4 2 



(6.37) 



(6.38) 



In this case we have 



o ldX A^) l0g (^ 



7 / 1 d ^ 72 2 lo S — [k -V 

t k l — y l s v 



The integral is awkward and may be made easier by using the following relation 

1 



~r~ log (A; + j/) = ^-Li 2 f^- (k — y) 
k — y dy \2k J k — y 



log (2fc) 



(6.39) 



(6.40) 



Substituting back in for fc 2 = | + \5 -f- J etc and dropping terms O (5 2 ) the integral is 

found to be 



Jo 1 da; A7r7 lo s 
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Li 2 (l + ^)-Li 2 



(6.41) 



This simple step is illustrated to condense the structure of divergences that are present in the 
limit m ej mf — > 0, the collinear divergences. A similar decomposition to that in equation 16.381 
gives the first two terms in 16.351 as 



Jo 



dx — ( log f — 

Ale 6 1 s 



2Tt[l- log [- s 



2 >g(-^) + log(-f 



(6.42) 
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Finally, taking 2k ~ 1 outside logarithms gives 




Su (t) = -TTTTT lim 

32ir z e'^o t \ \ t I ~ \ t I I \e ' \ s I l Jo £\{t) ' \ s 

(6.43) 

In summary we have that in the most general case the box diagrams under consideration have 
an amplitude which consists of a scalar IR divergent term 2ie 2 tSu (t) Asom (—2ie 2 uSu (u) Abotu 
for the crossed box diagrams) along with other infrared finite terms the exact nature of which 
depends on the physics of the exchanged X° boson. In addition it is possible that ultraviolet 
divergences are present among the other non-universal terms, in the case of the quantum gravity 
model of Giudice et al the amplitude is at most logarithmically divergent and we expect that 
divergence enhances such a term by a small factor. Adding the universal contributions of the 
four diagrams El together we have 

Abox | u niversal 

= Am + A D2 + Aci + A C 2 



lim f /_ 



H + ^ (£)) log (±) + I Jo 1 dx ^ log (Ml) - I jj dx ^ log 



A B orn 
(6.44) 

The integral final integral in 16.441 is the same as the one preceding it (shown in EES} with the 
replacement t — > u. These two integrals add to give terms which are finite as m e ,m/ — > 0. 
It is known in the literature that the initial state-final state interference corrections contain 
no collinear divergences, here we see explicitly that this result is true also for the universal 
correction as one would expect. We can safely take the limit m e ,m/ — > leaving 



Abox\u niversal ~ 



4 W < U I ,<W-) - i i<W± 
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(6.45) 



Denoting lim^o [••■] as JF, the interference of the box diagrams with the born amplitude may 
be written 

it , A_ /it \ _ 2a T V /it 



^spins [-^Born-^Box ^Box^-Born) — 3~ ^spins ■^-Born^Borni (6.46) 
\ 'IX 

hence 

2a 

d°Box\ Universal = —T&OBorn- (6.47) 

If we now add the differential cross section due to box diagrams interfering with the Born level 
amplitude to the differential cross section due to soft bremsstrahlung diagrams we obtain the 
total differential cross section for the universal part of initial state-final state interference as 

do~Universal = d(T Softjnt + Box \ Universal 

= V (log (^) ^g (f ) - Li 2 (l + f ) + Li 2 (l + i) . (6.48) 

+ >g 2 (^)-|l0g 2 (^))da B orn 

This is the universal contribution to the differential cross section from initial state-final state 
interference contributions to some new physics exchange process (X°) with tree level differential 
cross section dcrg orn . For s/s ~ 200 GeV while one may take the soft photon cut-off to be 
uj ~ 50 MeV (i.e. typical LEP experiments) this gives log (-^rj = 15, this constitutes a so- 
called large logarithm. This large logarithm log f^s) arises from the cancellation of the infrared 
divergences in the box diagram and initial and final state bremsstrahlung contributions. As 
noted earlier only the universal scalar integral terms are infrared divergent thus no other terms 
will undergo this large logarithmic enhancement. These universal factorizable corrections maybe 
resummed (exponentiated) as in the conventional treatment of the Sudakov effect. Considering 
only terms O (a log (jfjsj) corresponds to the leading log approximation for initial state-final 
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state interference. We stress that terms which are not universal are not infrared divergent 
and so they will not have the large logarithm log (4^2), they correspond to sub-leading terms 
O (a). We conclude that initial state-final state interference corrections are universal within the 
leading log approximation and we anticipate, conservatively, that the leading log (universal) 
terms dominate non-universal sub-leading terms by a factor of 10. 

It is of note that the theoretical study of the ISR*FSR interference by the author began with 
the calculation of the QED box diagram (i.e. for the case that X° is simply a photon). The 
FeynCalc [21] package was used to assist Dirac traces and reduction of the Passarino Veltman 
functions [22]. The universal leading log term was compared to the non-universal terms in this 
case and was found to dominate the sub-leading terms by a factor of at least 13 across the cos 6 
region under study. 

7 Summary and Conclusions 

We have seen that the radiator functions describing ISR, the dominant correction, are merely the 
result of an intrinsically factorized DGLAP type analysis with no additional non-factorizable 
components and that these describe the total cross section to better than the 1% level. In 
addition we have established that even ISR*FSR interference contributions to the differential 
cross section are universal to better than the leading log approximation (note that equation 
16.481 contains sub-leading terms which are nonetheless universal). With this in mind we feel 
justified in using the numerical radiator functions obtained from the ZFITTER program or its 
relatives to correct the tree-level differential cross sections of processes involving the exchange 
of some new particle. This method is valid within the context of the leading log approximation 
(and also somewhat beyond). 

We would like to stress however that though we conclude that to good approximation the 
radiative corrections we have discussed are independent of the new physics in the diagrams they 
are not independent of the topology of the diagrams. The ZFITTER package and this analysis is 
only strictly valid for the case of s-channel difermion processes. ZFITTER treats the t-channel 
processes i.e. e + e~ — > e + e~ differently to the other difermion final states. Our analysis of 
the DGLAP electron structure function approach to ISR will still hold for t-channel processes, 
these corrections pertain only to the incoming particles. On the other hand we expect that our 
analysis of the universality of ISR*FSR interference as discussed in section 16^21 will be modified, 
certainly one should expect the angular dependence of the radiative corrections to be different. 
Nevertheless a large logarithmic, universal, factorizable term will result. This term will be, 
by analogy to the s-channel result, of the form ^ log (4^2) log (...) i.e. the t-channel result 
should still involve a logarithm of the ratio of the hard scale physics scale to the soft scale (oj 2 ). 
Provided that the cuts on the cos 9 region are not too loose i.e. provided all the events in the 
sample are such that t 3> to 2 , the resulting enhancement for the universal term (log (7^2)) 




will still be large thus making the use of radiative corrections to e + e — > e + e in the standard 
model valid, to good approximation, for new physics e + e~ — > e + e~ processes. 
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Figure 7: The diagrams for tree level difermion processes. From left to right these are s- 
channel and t-channel processes respectively. Only electron-positron final states may result 
from t-channel processes. 

e + (p + ) f(qj e + (p_) f(qj 




e(pj f(q_) e(p + ) f(q_) 



Figure 8: Diagrams corresponding to difermion events involving bremsstrahlung. The black 
circles denote all possible one-particle-irreducible diagrams. The arrows denote the flow of 
fermion number, momenta p + , p_ are flowing inward, q+,q~ are flowing outward. 

Appendix A 

A difermion event is basically an event where an electron and positron in the initial state 
exchange a Z°* or a photon 7* and leaves a fermion and an anti-fermion in the final state. 
This can happen via s-channel t-channel exchange, these are depicted in figure at tree level, 
t-channel processes are only possible for a dielectron final state. 

Radiative corrections to these processes are significant in the experiment and the definition 
of a difermion process is modified. More generally we define a difermion event as an event 
where an electron and positron emit photons in the initial state then undergo an interaction 
which produces a fermion anti-fermion pair and photons, finally the fermion and anti-fermion 
may emit additional bremsstrahlung photons. Generally this can be reduced to three types 
of diagram, the tree level diagrams with bremsstrahlung from the external legs and vertices 
replaced by effective vertices (one-particle-irreducible diagrams with possible photons emitted 
from them) and one with photons radiated off the external legs and a single 4-point effective 
vertex, to lowest order this is a simple box diagram. In addition it is possible that photons are 
radiated off the one-particle-irreducible diagrams representing the vertices. 

We define the Mandelstam variables s, t , u as^ 10 ) 

s = (p++p_) 2 = (g+ + g_) 2 

t = (p--q-f = ( P+ -q + ? ■ (A.l) 
u = (p_ - q+f = (p+ - q-f 

Neglecting the mass of the electron it has four momentum p_ = (E b , 0, 0, E b ) prior to any 
( 10 ) Momenta pi are flowing inward, qi flow outward. 
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bremsstrahlung, this defines the z-axis for the experiment, the positron has four momentum 
P- = (Eb, 0, 0, — Ef,), consequently 

E b = (A-2) 

The initial state radiation (ISR) is defined as that which is emitted from the external legs of 
the diagrams. We denote the combined outgoing momentum of the radiation emitted from the 
external legs of momentum pi by kj (i — +,—, j = 1, 2) and from the external legs of momentum 
qi by kj (i — +,—, j = 3,4). Primed variables are defined to enable us to discuss the difermion 
process in the presence of radiative corrections where the invariant mass is reduced 

p'_ = p_ — k\ 

q~ = q~ + k 3 
q'+ = q+ + h 

In the presence of such radiative corrections there is some ambiguity in the definition of s, t, u 

as 

{p-+p + f ^(q-+q+) 2 (A.4) 
etc, hence for the case of radiative corrections we use 

s = (p_+p+) 2 
t = (p_-q_f 

u = (p+ - q- v 



,2 / \ 2 

/ I nrJ I nrJ \ I „l I „/ 



s = [p_+p + j = [q_ + q 

t' = (p'--q'-) = (V+ -<?:'" 

u> = (pL-q'Y = (p' + -q'y 



(A.5) 



Appendix B 

In this appendix we show briefly the integration of the soft photon phase space factor associated 
with the soft bremsstrahlung discussed in subsection 16.11 The integral to be performed is 

1 / 1 \ e ' r r\ £ ' +3 h 1 

(B.l) 



4ii \n ) J\k\<u) ko (k.pi) (k.pj) 
Following 't Hooft and Veltman (2jj we introduce a parameter k and define 



V = KPi , (B.2) 
q = Pj 

where k is defined to be the solution of (p — q) 2 = which gives go with the same sign as po — qo- 
With these definitions lB.il becomes. 

-(-Y f ^ 1 (B.3) 
4ir \fi J J\k\<Lj k (k.p) (k.q) 

now combine the denominators with a Feynman parameter (afe) -1 = Jq dx [ax + 6(1 — %)]~ 2 - 
This easily gives, 

4tt \n) J \k\<uj k JO ux (Epf 

4*\p) JO aj -J\k\<u k k 2( Po _\p\ cos gy 
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where we have defined P = q + x (p — q) . 

The integration measure can be written explicitly in terms of its angular components and 
momentum component in the usual way. The angular integrations are simplified by essentially 
redefining the photon momentum space axes such that 9 is the polar angle, in this case all but the 
polar angle integration is trivial. This is analogous to how integration over the azimuthal angle 
in three dimensional problems gives a factor 2ir when the integrand has rotational symmetry 
about one axis, the d-dimensional analogy is a well known result see e.g. [S] 
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(B.5) 



Equation IB .21 gives the integration measure dO^ on a o?-sphere, in our case the phase space is a 
d = n — 2 dimensional sphere of radius k embedded in n — 1 dimensions. Hence we decompose 
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-d cos 9 this becomes, 
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Expanding in e', dropping terms O (e') and above this becomes, 



K 



2 Jo 



dx 



+i 



dC 



Po 



P 



C 



1 n w , 7 
- + log - 
e' ix 



_^_l 0g7r+ _l 0g (l_C 2 



(B.7) 



(B.8) 



with cos 9 relabeled as C. The angular integrations maybe safely carried out using a symbolic 
computer algebra package giving for integral IB.8I 
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(B.9) 



(recall P = q + x (p — q)) where we have introduced the infrared regulator 

1 2 

- = —, + 7 + log 7T. 

e e 

Again by analogy to [22], the integration is transformed such that it is over m = Pq 
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where we have used the fact that k was defined so that (p — q) = in determining P and 

- 2 

P above. In terms of the new variables we have 
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The rest of the integrand and the integration measure transforms as 
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Differentiating lB~T2l with respect to m we find 
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this enables us to rewrite IB. 131 as 
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P 2 



P 



log 



Po- 


— * 

p 


-]- 2 


Po + 


p 


/ m — 
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2vl \ dm 

Substituting these transformed qn a nt it ies IB . 1 31 and IB . 1 51 into the phase space integral IB. 71 gives 

d £, + 3 fc l 

<u k (k.pi)(k.pj) 

(B.16) 

log 



± 1 

4-7T \/l 
K ( 1 

2 \ e 



log (a) 3 + 2 log 2) J? dx^, 

I K [P0~\p\ dm ( d l 
"'" 2 J<20-|?1 2uZ I dm U ' 



V ' m— i(2ii(jo— <3 2 ) y m— d 

The term log (...) log (...) can be rewritten as |t- log 2 
trivial. The first integral is also easy, 



/ X i^ = 2^ 



2^ 

g 2 



in which case the integration is 

(B.17) 



This leaves one integral 



d e '+ 3 fe 1 



Lit fc <w fco (k.pi)(k.pj) 

2 + 21og2) l g(l + ^) 



+ 



2d 



m(m—v) 
q 2 +mv—2vqo 



(B.18) 



f/0 _ k1 2,/ J»-|5| am 



log 



m—v m 

v ' m-£(2i;go-g 2 ) 



The final integral is awkward and is done with the help of the following two identities 



log (m — a) 
£Li 2 (A (m - „)) 



log (l - ^) + log (t; - a) 
3— log (1 — A (m — u)) 



rPo-b1 dm _^io g ( 2iz2._ 



-Li 2 (1 - f ) + Li 2 ( (4^ 



Po-|p| 
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+ 



rP0-\p\ - 1 i f (m-v)v 



(B.19) 



(B.20) 



Finally we abbreviate v (v 2 — 2vq + q 2 ) 1 
x = (m — v) c, dm = -dx hence 



c and perform a transformation of variables 



v c 



dx ± log (x) 

1 l 2 ( i>(m— ■») 



PO-\p\ 

qo-\q\ 



(B.21) 



For the phase space integral f6. Ill we now have 

d s '+ 3 fc i 



k \<U! k (k.pi)(k.pj) 
^(l+ 1 °g(^)) 1 °g( 1 + ¥) 



+ 



4vl 

K 

2vl 



>g 2 



m(m— v) 
q 2 +mv—2vqo 



v(m—v) 
q 2 +v 2 —2vqo 



+ Li 2 1 



v(m— v) 
2vqo,—q 2 —v 2 



Using the identity for the dilogarithm identity Li 2 
becomes 

-k J_ flV' r d E '+ 3 fc l 

47T U, 



Po-\p\ 

qo-\q\ 
(B.22) 

-Li 2 (z) - Mog 2 (-z) -C(2) this 



|fc|<oj feo (k.pi)(k.pj) 

^(! + log(^))log(l + ^) 



(B.23) 



+ — 

~ 2uZ 



i log 2 (; 



m(m—v) 



Li* 1 - a 



Li 2 f 2 *W-" 2 



PO-|p| 
90— 1^1 



Retracing the algebra of IB .121 backward and using the relation p 2 — 2wp = q 2 — %vq we can 
write 



(go— |g|)(go— |g|— 

q 2 + (qo-\q])v-2vqo 

(po— Ip1)(po— IpI— 
p 2 +(pO-|p|)«-2'Upo 
2vqg— q 2 — v 2 

2vpo—p 2 —v 
v(po-\p\-v) 
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90+ 
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PI 


P0 + 


Pi 



^— go— |g| 
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^-po-|p] 
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Also for any momentum four vector V, in the limit of small masses, V 2 <C Vq and we can 
approximate, 



V 
V 



Vo-V (l 



V + Vn 1- 



V 1 



2V 



(B.25) 



Specifically in the case of initial state-final state interference contributions to dasoft we re ~ 
quire four phase space integrals corresponding to all possible permutations of a bremsstrahlung 
photon from an initial state leg and a bremsstrahlung photon from a final state leg i.e. all 
combinations of pi and pj where pi — p_, p + and pj — q~, q+. Working also in the limit that 
the photon carries off no energy p^o = Pj t o, the phase space integral can now be written 



I (J) f\> 



d € '+ 3 k 1 



fc|<u; feo (k.pi)(k.pj) 

ti{i + 1^ (*?)) log (l + f) + £i log 2 ) - & log 2 

^— mo— N 



Itt 

K 
K 

2vl 



Lio 



u— mo + | - m| 
v 



+ Li 2 



m=Kpi 
m=pj 



(B.26) 



Returning to equation ESI we see that the soft photon contribution to ISR*FSR bremsstrahlung 
is 



da 



a 



Softjnt 



d 3 k 



t 



+ 



t 



u 



4n 2 J\k\ <L0 k \(k.p + ) (k.q+) (k.p-) (k.q-) (k.p+) (k.q-) (k.p-) (k.q+) 



da Born 

(B.27) 
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Finally we must determine k for each phase space integral where n was defined earlier as the 
solution to 

{m - Vif = 0, (B.28) 

for which Kpi$ —pj,o and p^o have the same sign. We shall demonstrate how to obtain k for the 
case pi = p + pj = q + , the results generalize easily to the other combinations of momenta. In 
this case expand ing IB . 281 gives a simple quadratic equation for k. Working in the limit t 3> 
gives the two solutions ( ) for k to first order as 



t 



(B.29) 



In the massless limit the Mandelstam variable t is 



t = -|(l-cos0++), (B.30) 

where 9 ++ is the angle between and \q+\, therefore t is a negative quantity. We require the 
solution for which sign (kp +>0 — q+,o) = sign (g+,o) *-e. 

sign Qv 7 ^ (k - 1)) = sign C^V^j (B.31) 

so we require k > 1. Again working in the limit t S> m 2 f this means taking solution k + = — ^ . 
Exactly the same mathematics and the same k are obtained using pi = P-, Pj = q~ and only 
marginal differences in working give k + = — \ for the other two combinations of momenta. 
Noting that k is large we have that to lowest order in small things 



= 1 «2 j? fa _ _g \ fa _ JUL 



(B.32) 



Now considering solely the dilogarithm terms in the small mass approximation IB.25[ we have 
for k = — K 



Li 2 fl- m °- H +Li 2 l- m ° + H 



l m= * 1 (ml\ / s\ vr 2 



V 



m=pj 

> • 



■log - -L H 1 , 



(B.33) 

where we have used Li 2 (-^) = -Li 2 (-^) - § log (^) -£ (2) and Li 2 (1) = £(2). The same 
term for the second solution k = — A- takes the same form as above but with the replacement 



t — > u. Substituting all of this into IB .271 gives finally 

d. s „ /M „« = v [(1 + log (¥)) log (s) " Li2 (' + f) + Li2 ( : + s 
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